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ABSTRACT 


An experimental investigation has been carried out on an underexpanded jet of air 
issuing from a converging rectangular nozzle of moderate aspect ratio. Schlieren pictures 
of the flow field along with hot-wire measurements in the jet were obtained at different 
pressure ratios. At the pressure ratio corresponding to the maximum screeching sound, 
Schlieren photographs show a very strong organized cylindrical wave pattern on either side 
of the jet, with their respective sources being located at the end of the third shock cell. 
Associated with this wave pattern is a large increase in the angle of spread of the jet. It is 
shown that the self excitation helps to induce large-scale vortical motions in the jet both 
in the plane containing the small dimension of the nozzle and in the plane containing the 
long dimension of the nozzle. However, the locations of these structures are different in 
the two planes. Nevertheless, the characteristic Strouhal number corresponding to these 
large-scale structures in both planes is the same and equal to 0.12. The influence of the 
self excitation on the mean velocities and rms intensities has also been investigated. 

For the full range of pressure ratios studied, similarity was found both in the mean 
velocity and rms intensity profiles in the two central planes beyond 80 widths downstream 
of the nozzle exit. However, the shapes of the similarity profiles are different in the two 
planes. For the downstream distances studied, complete axisymmetry (identical mean 
velocity profiles in both planes) was not found, which suggests it may persist for a large 
distance downstream of the nozzle exit. 
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INTRODUCTION 


The increased interest in V/STOL aircraft has caused a great deal of attention to 
be/focussed on jet flows exiting from rectangular nozzles. In particular, the mixing 
characteristics of an underexpanded jet exiting from a converging rectangular nozzle of 

moderate aspect ratio is of great interest. As part of an ongoing program at Stanford to 

study the mixing characteristics of single and multiple rectangular jets both in free and 
confined configurations, a systematic investigation has been carried out on the structure of 
a jet issuing from a moderate aspect ratio rectangular nozzle at flow speeds ranging from 
low subsonic to supersonic. The detailed measurements of the mean and turbulent velocity 
components of an incompressible jet was reported by Krothapalli, et al., 1081. The effect 
of the exit Mach number on the overall flow field of a subsonic compressible jet was studied 
by Ilsia et. al., 1082a. This paper reports the results of a study of a supersonic rectangular 
jet (choked at the nozzle exit). In all these investigations the nozzle used was the same, 
thus providing an opportunity to assess systematically the various changes brought about 
by the varying flow conditions. ... ... ... 

It is well known that the structure of a choked undcrexpanded jet displays different 
features to those of subsonic and ideally expanded or shock free supersonic jets. These 

features include discrete tones in the sound spectrum, known as screech tones, under certain 

conditions and an accelerated spreading of the jet with enhanced mixing. The jet screech 

phenomenon, an oscillatory condition common to underexpanded jet flows, is described by 

Powell 1953, as arising through a feed-back mechanism. A disturbance in the jet shear 
j layer is convected downstream (assumed amplified in passing downstream) and strikes a 

: cell boundary, scattering intense sound at that point. This sound propagates through the 

• ambient medium in the upstream direction and interacts with the jet shear layer close 

r to the jet exit, giving rise to a new downstream travelling disturbance that continues 

I the cycle. The essential element of the cycle is the amplification of the downstream 

f travelling disturbance. This process maintains itself mainly because the energy lost through 

| .* scattered radiation and viscous effects is resupplied to the propagating disturbance by an 

I instability of the flow. The detailed mechanism of screech tone generation is yet to be 

I . understood, although the general nature of the feedback loop, was treated by Powell who 
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developed a simple relation between tone frequency, cell size, and operating conditions. 
Under certain conditions, the concomitant sound production accelerates the spreading of 
the jet and enhances mixing. However, the conditions for the existence or absence of such 
a phenomenon cannot yet be predicted. 

A considerable amount of experimental data is available in the literature to charac- 
terize the near sound field of a choked undorexpanded jet (Powell 1053a, 1053b; Lassiter 
and Hubbard 105G; Davies and Oldfield 1062; Hammilt 1061; Harper-Bourne and Fisher 
1073; Westley and Woolley 1073; Sherman, Glass and Dullep 1076; Krothapalli, BaganofT, 
Hsia and Karamchcti 1081, 1082). Most of the work reported in the lit erature has been as- 
sociated with the flow field of an axisymmetric jet. Although the structure of a rectangular 
jet exhibits some features similar to that of an axisymmetric jet, there exists some impor- 
tant differences between the two flows. Some of these differences even make experimental 
study easier. For example, the nearly two-dimensional structure of the disturbances and 
their associated sound fields near the nozzle exit are easily identified in flow visualization 
pictures because of the long optical path length present. It is thus experimentally attractive 
to study the interaction between the sound field arid the mean and turbulent components 
of the flow for the flow exiting from a rectangular nozzle. 

The effect of acoustic feed back on the spreading and decaj of an axisymmetric 
underexpanded jet was first studied by Glass 1068. However, his experimental study was 
limited to the measurement of gross properties of the flow because of his use of a pitot tube » 
in the experiments. None of the works that followed presented information regarding the 
effects of screech tones or organized sound waves on the turbulent structure of the flow; in 
particular, none on the rectangular jet, the objective of the present paper. The detailed 
structure of the near sound field and its source for an underexpanded jet operating at 
maximum screech condition is to be given in a separate paper by Krothapalli ct al., 1982. 

The principal parameters or variables governing the flow of a free rectangular jet 
are the pressure ratio R (stagnation pressure/ambient pressure), Mach number A/, and> 
Reynolds number Re, the turbulence level of the flow at the exit of the nozzle, the 
conditions of the medium into which the jet is issuing, and the aspect ratio of the nozzle. 
In the present investigation the pressure ratio, R, was varied from 2 to 5.8. This interval 
corresponds to a Mach number range, based on fully expanded isentropic flow, of 1.05 to 
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1.8. The Reynolds number employed here is based on the width D of the nozzle (small 
dimension) and given by Re — McD/u, where a and v are the speed of sound and 
kinematic viscosity of the ambient medium respectively. This Reynolds number was varied 
from 7.2 X 10 4 to 1.0 X 10 s in the experiment. A rectangular nozzle of aspect ratio 10.7 
was studied. The inlet geometry' of the nozzle was designed specifically to obtain a low 
turbulence level at the exit plane of the nozzle. The total-pressure profile at the exit plane 
of the nozzle was found to be quite flat. 


3 



APPARATUS, INSTRUMENTATION AND PROCEDURES 

A high-pressure blow-down type air supply system was used to provide the air flow to 
a cylindrical settling chamber having dimensions of 1.75m long and 0.6m in diameter. The 
temperature in the settling chamber was maintained at a constant temperature, usually at 
room temperature, to an accuracy of about 0.5°C over the duration of each test. Before 
reaching the nozzle, the air was passed through an adapter, containing six screens set 
5cm apart, to minimize disturbances at the nozzle inlet. The ratio . reas between the 
adapter and the nozzle exit was about 40. Two nozzle sizes were usc^. .a the study. The 
dimensions of the rectangular exit of the small nozzle used was 50mm long (L) by 3ir*n 
wide ( D ). The exit dimensions of the large nozzle were 50mm long and 5mm wide. The 
nozzle exit in each case was preceded by a 40mm long smooth rectangular channel (50mm 
X 3mm and 50mm X 5mm respectively). The 3mm nozzle used in the investigation, shown 
schematically in Figure I, was a single central lobe of a multi-Iobed nozzle employed in a 
related investigation. Most of the detailed measurements were made using the 3mm nozzle. 

The experiment made use of the same model and air supply system as that described by 
Ivrothapalli et ah, 1981. A 7mm thick acoustic insulation was used to cover all exposed 
surfaces on the nozzle to minimize acoustic reflections from these surfaces. 

A conventional Schlieren system was used for flow v.sualization purposes. The configuration 
employed was a single pass design, with the optical axis folded twice using two 25cm 
diameter, 3.05m focal length spherical mirrors. The light source employed was a strobo- 
scopic flash unit having a 1.5/r sec flash duration. The stroboscopic feature was used for 
visual observations while single pulses were used for photographic records. 

For downstream distances X/D less than 15, where the velocities are high and the 
flow contains shock waves, measurements were confined to the use of a pitot tube having 
an external diameter of about 0.40 mm and a diameter of the sensing tube of about 0.25 
mm. Further downstream, in regions where the flow is known to be subsonic, hot-wire 
and/or hot-film anemometry was used. 

All velocity measurements were made with a DISA 55 M10 constant - temperature 
anemometer in conjunction with a DISA 55 M25 linearizer. The hot wire used was a 
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DISA 55 Pll single normal platinum-coated tungsten wire with a f/t m diameter and 
lmm length. The hot film used was a DISA 55 R31, a wedge shaoed 0.5/4 m quartz 
coated nickel sensor with a 0.2mm width and 1mm length. The frequency response of 
the hot-wire electronics, in response to a square wave test with the probe placed in a 
high speed (Af = 0.8) low turbulence jet, was approximately 40 KHz. The hot-wire and 
hot-film probes were calibrated using a low turbulence air jet regulated by a DISA 55 
D0O calibration unit. The flow velocity, with correction for compressibility, is obtained 
as an analog voltage. The hot-wire output voltage was linearized for the velocity range 
20m/sec to 280m/sec, with the output voltage E directly proportional to the velocity 
U{E = KU , where K is a constant selected in conjunction with the linearizer gain control). 
Similarly, the hot film output voltage was linearized for the velocity range 20m/sec to 
300m/scc. Despite the direct proportionality between the fluid velocity and the hot-wire 
voltage during calibration, the hot-wire fluctuation voltage r, in a turbulent compressible 
flow, will depend additionally upon fluctuations in the density and the total temperature 
(Horst inann and Rose, 1977). Lacking proper means to accurately correct velocity data for 
these fluctuations, rms values will be expressed in terms of e =< c > /E (where < e > 
denotes rms), which becomes equal to the turbulence intensity < u > fU for low Mach 
numbers. The detailed calibration procedure used is given by Hsia et al., 1982b. 

The signal from the linearizer was passed through a DISA type55D3l digital voltmeter, 
a DISA type 55D35 rms unit, and a TSI model I07G voltmeter to get the mean and rms 
values. Integration times with these instruments can be selected in discrete steps from 
0.1 to 100 sec, and 10 sec was typically used. Spectral measurements were made using a 
Nicolet type GG0B spectrum analyzer. 

A cartesian coordinate system (A',V, Z), defined in Figure 1, was employed with its 
origin located on the centerline of the jet. Hot-wire or hot-film traverses were made in 
the two central X, Y and X,Z planes at streamwise locations (A') from 40 1) to IG0D. In 
all the measurements, the sensor wire was oriented parallel to the Z axis. Mean velocity 
measurements were made across the entire jet in order to establish the symmetry of the 
flow about the central planes, however, only the data for each half plane will be presented.' 

The controlling parameter in this investigation was the stagnation pressure, p 0 , which 
varied from 30psia to 85psia, and was maintained within an accuracy of ± 0.02psia. This 
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range corresponds to nozzle pressure ratios (stagnation pressure/ ambient pressure) varying 
from 2 to 5.8. The temperature in the settling chamber was maintained to within ±0.5°C 
of the ambient temperature. Experiments were conducted for a number of pressure ratios 
using both nozzles, however, only a limited selection is presented here. The detailed hot- 
wire measurements were made using the 3mm nozzle and were limited to the three pressure 
ratios R = 2.7, R = 3.8, and R = 5.4. 



RESULTS AND DISCUSSION 

(a) Flow Visualization Stud) 

/ ■ ‘ ' . 

Typical Schlieren pictures of the flow issuing from the 3mm nozzle for three different 

pressure ratios are shown in Figure 2. The knife edge was oriented parallel to the jet 
axis in these photographs. The Schlieren system was purposely adjusted to enhance the 
view of the wave structure in the outer region, rather than to shew the details of the jet 
flow. Consequently, features such as shock cell structure are not seen clearly in these 
photographs. The long optical path associated with the rectangular geometry make- 
possible t' ' viewing of wave systems that otherwise would not easily be seen These 
pictures display the flow field from the nozzle exit to a down-stream location of about 

30 D. 

Pqj* pressure ratios i? ^ 1.0, the flow beyond the nozzle exit is supersonic and results m 
the formation of a series of shock cells, as shown in Figure 2. The most striking observation 
for the three cases shown’ in the figure are the large changes in the angle of spread of 
the jet, the appearance of a very strong organized wave pattern, and the development 
of a M <eh wave radiation pattern at the highest pressure ratio, /or example, the total 
angle of spread of the jet is about 20° for both the lowest (It = 2.7) and the highest 
(ft = 5.4) pressure ratios, while the total angle of spread is about 30° at the pressure ratio 
for maximum screech sound radiation (/? = 3.8). Associated with the large spreading 
rate is an organized cylindrical double wave pattern which originates alternatively from 
the two sides of the jet. The center of the arcs formed by the sound waves on each side 
can be located and identified as the source for each wave system. This source is loca.cd 
approximately at the end of the third shock cell (x/D = 7). The picture seems to suggest 
that a principal wave length exists and that this wave length '.s the same for all directions 
of propagation of the wave. The second wave front seen in. the figure introduces a second 
and much shorter characteristic length, which corresponds to a frequenc> much too high 
to record with ordinary detectors. A similar double wave pattern was observed by Powell 
1053c. in his edge tone excitation of rectangular jets. For the pressure ratio R = 5.4, 
high frequency sound waves, emanating at an angle of about 45° to the direction of the 
stream, can be seen radiating from a region close to the jet exit. This apparent, radiation 
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from flip supersonic region of the jet was nlso observed by bowson nnd Ollerheiul IflfiS, 

• ’ - * . 

in their experiments with ideally expanded axisymmetric jets, and they identified it with 
Mach-wave radiation. From the near, field sound spectra, not shown, here, discrete tones, 
known ns screech tones, are observed in all three enses (Krothapnlli et nl., 1081). However, 
the simultaneous existence of an organized wave pat tern with large spreading appears only 
in the range of pressure ratios from 3 to 4. On reviewing all the pictures taken in the 
study, an organized vortical structure appears to he identifiable in the shear layers on 
either side of the jet, and these grow in intensity as the pressure ratio approaches the 
value for maximum screech sound generation, which occurs here at It — 3.8. 

To study the degree of two-dimensionality present in the How, j * otographs were taken 
at right angles to the view shown in Figure 2. These results are presented in Figure 3 
for the same conditions as those in Figure 2. Three shock cells can be identified in the 
photographs by the three bright lines in the jet, which are the shock waves at the end of 
each of the corresponding shock cells. For downstream distances r less than It)/), organized 
lines are observed lying along the span of t ho jet . The spacing between them seems to 
increase with downstream distance. These lines suggest the presence of spanwise vortieies 
developing in the shear layers surrounding the cell structure (see Figure 2). Uecause of 
the concave nature of the shear layers, it is suggested by Morkovin 1081 that they may be 
Cloertler vortieies (Schlichtiug 1080). These features are more clearly seen in Figure 4 to 
be. discussed below. Wry little spreading is observed for the lowest and highest pressure 
ratios, while at It = 3.8 the jet width first reduces at intermediate distances and then 
increases with downstream distance. 

Figure 1 presents an example of the flow from the Jimm nozzle for conditions close 
to those for Figures 2b and 3b. In this case the aspect ratio of the nozzle is lower (10) 
than for the 3mm nozzle (10.7). On comparing the two flows in the plane containing the 
small dimension of the nozzle, one finds that a single wave pattern is present in Figure 
la, while a double wave pattern is present in Figure 2b. This feature was found in all of 
our tests and appeared to be associated with thy larger aspect ratio of the 3mm nozzle. 
'Phis may be consistent with the fact that when the aspect' ratio is further reduced no 
double-wave pattern is seen, as in the case of circular jets (Thomson 1071). Uecause of 
the larger nozzle size, the shock cells can be seen more clearly in Figure 4a and evidence 
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of their large amplitude lateral oscillations is quite apparent in the photograph. Likewise, 
the source of the strong cylindrical waves seems to be located at the point where the third 
shock cell breaks up. The large scale vortical structures mentioned above can be seen more 
clearly in the picture and can be followed by their orderly dark patterns. 

’ ■ . . ' i 

- - --• . t • • 

The photograph taken at right-angles to the view in Figure da is shown in Figure 4b. 
Three shock cells can be identified in the figure, and very little spreading of the jet occurs 
in this view, indicating the flow is fairly two dimensional within the region covered in 
the photograph. The organized spanwise disturbances or vorticies can be seen even more 
clearly in this photograph and can be followed for a considerable distance downstream. 

It has been shown by Poldervaart et al. 1973, that the placement of a reflective 
surface near the exit of an underexpanded rectangular jet has an important effect on the 
flow. To further investigate this effect, several reflector geometries were used. Two typical 
geometries are shown in Figures 5 and 6. The conditions were the same as those in Figure 
4. The most striking conclusion drawn from the experiments is that the flow can be either 
destabilized, characterized by a large spreading of the jet, or stabilized, by use of a small 
plate, and the effect is present with the plate located on only one side and in almost any 
position. Figure 5 gives such a comparison where the distance between the plate and the jet 
was adjusted while keeping the plate parallel to the jet axis. In this case, the distance the 
plate must be moved is essentially related to the wave length of the sound wave produced. 
Although similar stabilizing effects have been observed by others, e.g. Hammitt 1961. with 
acoustically absorbing baffles at 90°, we believe this is the first case observed where the 
effect is so pronounced when using a hard reflecting surface. Figure 6 is an example where 
the plate was positioned 135° to the jet axis. The presence of the plate seems to generate 
an additional wave system whose origin is located at about the end of the first shock cell. 

(b) Total Pressure Measurements 

To make measurements in a compressible flow with both a high turbulence level and 
the presence of shock waves is difficult. This being the case in the near region (t/D < 15) 
of the jet, a pitot tube was used to measure the total pressure and establish some of the 
basic features of the jet. 

A typical variation of the pitot pressure, 7 ) 02 , along the centerline of the jet for a 
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pressure ratio of 3.8 is shown in Figure 7. The absolute values of po2 are divided by the 
ambient pressure P a and plotted against the nondimensionalized downstream distance. 
These results cannot be converted directly into velocity or Mach number because of the 
unknown entropy increase that has occurred. In steady supersonic flow with a single 
normal shock wave ahead of the pitot tube, a large pitot pressure correspond to low Mach 
number, and vice versa. Also shown in the figure is the location of the source for the 
sound waves observed in the Schlieren photograph (Figure 2b). Two shock cells are clearly 
identified by the cyclic behavior of the pitot pressure for x/D less than 6. The presence 
of a complete third cell is not evident here because of the time averaged data taken with 
the pitot tube. Phase locked Schlieren pictures taken within a single cycle (Krothapalli 
et al. 1982) show that a complete third shock cell is present only for a small fraction of 
the time during a cycle. It is observed that the downstream decay of the pitot pressure is 
considerably accelerated after about 8D, and subsonic flow is found beyond x/D = 15. 

The conditions at the nozzle exit play a controlling role in the downstream develop- 
ment of the flow, and thus these conditions should be specified in each experiment. In view 
of this, pitot pressure profiles in the central x,y plane for two different downstream loca- 
tions close to the nozzle exit (x/D — 0.1, x//; = 1) have been made. Figure 8 shows the 
absolute pitot pressure plotted against the position y. Close to the nozzle exit (i.e. for x — 
0.1/;) the pressure is nearly constant for two-thirds of the nozzle width. A typical profile 
within the same shock cell at x/D = 1.0 is also shown in the figure. If the pitot pressure is 
converted into Mach number, the profile (x/D =1.0) has a shape quite similar to the one 
observed within the first few widths of a planar jet. In many practical applications such as 
thrust augmenting ejectors for short take-off and landing aircraft and certain jet engines, 
the desired condition at the nozzle exit is a top-hat mean velocity or pitot pressure profile. 
In this experiment such a profile is obtained. 

(c) Mean Velocity Field 

While many quantities are of interest in identifying the structure of a turbulent jet, 
such as the turbulent stress tensor, various velocity correlations, dissipation terms, etc., 
some measure of the spreading rate of the jet along with mean velocity profiles are very 
useful for comparison with the behavior of other flows. Most of the mean velocity profiles 
were obtained using a hot wire in the downstream flow and a hot film in the regions close to 
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tin* nozzle exit. Both instruments gave consistent mean velocities (accurate to about 593) 
when used to make the same measurement. Measurements were made across the entire jet 
in some eases in order to establish the symmetry of the flow about its central plane. 

Figure 0 shows the measured decay of the mean velocity along the center line of the 
jet for three different pressure ratios, where U t is the measured center line velocity and 
f'o is the calculated velocity at the center of the nozzle exit (assuming ideally expanded 
isentropictropie flow). The data presented in the figure are for down-stream stations 
ranging from 40 D to 160Z). In all three cases for x/D greater than about 80, the center-line 
velocity decays as x -1 which is depicted by the solid line. This region is usually referred 
to as the axisymmetric region of a rectangular jet, where the velocity decays at nearly the 
same rate as that of an axisymmetric jet. Ahead of this region is a two-dimensional type 
region where the velocity decays at a rate roughly the same as that of a planar jet. 

The density of a jet relative to the density of. its surroundings has a significant effect 
on the rate of jet decay (Abramovich 1963). In particular, the greater the ratio of jet 
density to ambient density, the less rapid will be the jet decay. In the present experiments, 
the stagnation temperature was equal to the ambient air temperature and therefore the 
static temperature at the nozzle exit was less than ambient. For example, picking the case 
/? = 5.4, the ratio of the exit temperature for an ideally expanded gas to the ambient 
air temperature was about 0.62, and the ratio of the corresponding densities was therefore 
1.61. The effect of increased jet density is evident from a comparison of the data for the 
lowest and highest pressure ratios. However, for the case of R equal to 3.8 and for x greater 
than 80 D a more rapid decay of the velocity is noticed as compared to the lowest pressure 
ratio. This is principally the result of the strong acoustic excitation present downstream 
of the nozzle exit (see Figures 2b and 13) which causes an accelerated spreading of the jet. 

Figure 10 shows the distribution of mean velocity across the narrow dimension of 
the jet. in the x,y plane at different downstream stations, ranging from 40 D to 140/), for 
the lowest and highest pressure ratios tested. The mean velocity U is normalized at each 
station with respect to U e , and the distance Y is normalized by the distance (x — to) to 
the station in question, where xo is the distance between the virtual origin of the flow and 
the nozzle exit. For the case of R equal to 2.7 and for A' greater than 40 D, the profiles arc 
geometrically similar within the limits of error for the experiment. Self similarity is also 
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observed in the profiles for the case of R equal to 5.4, which is represented by the solid 
line (average of the data points not shown here). In the absence of any strong acoustic 
interaction with the jet, it appears that the self similar profile broadens with decreasing 
pressure ratio, i.c the spreading rate increases. From figures 0 and 10 it is observed that a 
wider profile corresponds to a smaller jet decay. For comparison purposes, the self similar 
profile of an incompressible rectangular jet of aspect ratio 38 (Gutmark and jgnanski, 
1076) is also shown in the figure and seems to agree well with the case R = 2.7. However, 
the presence of a strong organized wave structure in the near acoustic field of the jet, which 
occurs at R equal to 3.8 in this experiment, has a considerable effect on the distributions 
of mean velocity across the jet in thex.y plane. Such profiles are shown in Figure 11 for 
different down-stream stations ranging from 40 to 140 widths. As before, the distance y is 
normalized with respect to (x - x 0 ), where the value of x is different for the two segments 
40 < x/D < 80 and 100 < x/D < 160. As discussed below, the halfwidth (the width 
corresponding to the point where the mean velocity is equal to one-half of its value on 
the axis) in the x,y plane varies linearly with x, with two different slopes, thus producing 
two different virtual origins (see Figure 15). It is observed that the profiles in each of the 
two regions are geometrically similar except the scale factor is different for the two. The 
similarity profile for the downstream flow (.V > 100 D) is quite similar to that observed in 
Figure 10 for the lowest pressure ratio It — 2.7. The unusually rapid spreading of the jet 
seen in Figure.(2b) is the cause of the wider self similar profile in the region x < 80 1). 

Normalized mean velocity profiles in the x,c plane at different dow’nstrcam locations 
ranging from A0D to 160 D, and for the lowest and highest pressure ratios, arc shown in 
Figure 12. The distance z is again normalized with respect to the distance (x-xo) where 
x 0 is the location of the virtual origin (found by extrapolating linearly the variation of the 
half width in the x, z plane (see Figures 14 and 15). For the case of R — 2.7 and for r/I) 
greater than 80, the profiles are again geometrically similar, and have a shape quite similar 
to that of an axisymmctric jet. A corresponding observation is made for the pressure ratio 
R = 5.4 and is indicated in the figure by the solid line which was drawn through the data 
points (not shown here). It is observed that the similarity profile for the highest pressure 
ratio is narrower than that for the lowest. In light of Figures 10 and 12 it appears that 
the jot spreads at a slower rate with increasing /?, except in the range where self acoustic 
excitation plays a large role. This is consistent with calculations made for the ease of 
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a planar jet by Abramovich 1963. Figure 13 shows the corresponding normalized mean 
velocity profiles in in the x, z plane for the case of maximum sound radiation, i.e. R — 
3.8. As before, it appears that the profiles are geometrically similar for large downstream 
distances. The similarity profile in this case is wider than those shown in Figure 12. 

The growth of the jet in x ,y and x, z planes is shown in Figures 14 for the lowest 
and highest pressure ratios, while Figure 15 presents the case for maximum screech sound 
radiation. The ordinates j/jy 2 and Z|/ 2 are the respective distances from the centerline 
of the jet to the point where the mean velocity in each plane is equal to one-half of its 
centerline value. The jet in the x t y plane spreads linearly with x -and the locus of the half 
. velocity points is given by 

y,/ 2 = k{x - x 0 ); 

where k and x vary with pressure ratio. For the lowest pressure ratio the values of k and 
x are 0.105 and —5.3 D respectively. The corresponding val ;es for the highest pressure 
ratio are 0.08 and — 21.3£>. In the absence of any dominant discrete sound radiation, a 
decrease in the spreading rale of about 25% is found for an increase in the pressure ratio - 
from R = 2.7 to R — 5.4. We may note that for a subsonic rectangular jet, the constant 
k has a value of about 0.1, and does not change with increasing Mach number (Hsia et al. 
1982a). Similar observations were made by Maydew and Reed 1963 in their, investigation 
of turbulent mixing in axisymmetric compressible jets. 

It is known that an increase in the Mach number of a supersonic jet results in a 
decrease in the spreading angle of the mixing region at the boundary of the initial port ion 
of the jet. In most such experiments, the Mach number is increased by increasing the 
reservoir pressure which causes an increased flow density (if stagnation temperature is kept 
constant) and the observed effects were attributed to this increased density ratio between 
the jet and the external gas (Abramovich 1963, Glass 1968). A systematic investigation 
of the density effects in turbulent mixing layers was carried out by Brown and Roshko 1 
1974. Even though this study deals mainly with two dimensional shear layers, some of the 
implications of their results apply to the present investigation. 

The variation of the half-width in the x , z plane, shown in Figure 14, exhibits a 
linear increase for downstream distances x greater than about 80Z>. However, for x/D 
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less than 80 it shows a decrease. At some intermediate location the half widths in the 
two central planes 'cross each other. The distance from the nozzle exit to the crossing 
point increases with increasing pressure ratio, except for the special condition of maximum 
screech sound radiation. Hence, it may be suggested that the two dimensionality of a 
rectangular supersonic jet is preserved for longer downstream distances with increasing 
pressure ratio. 

Figure 15 shows the variation of the half widths in two central planes with downstream 
distance for the condition of maximum screech sound radiation ( R = 3.8). Except for the 
accelerated growth of yjy 2 for x/D less than 80, the overall behaviors of the half widths in 
both planes are quite similar to those observed in Figure 14. In the region 40 < x/D < 80, 
the constant k for the yj/ 2 variable assumes a value of about 0.16, while for x/D > 80, it 
has a value of about 0.1. Corresponding values of xq were determined in the two regions 
and used in the normalization of y (i.e. y/x — xo) in plotting the velocity profiles. It is for 
this reason that the velocity distribution in the x, y plane shows two different self similar 
profiles, as seen in Figure 11. The reason for the large spreading rate for x/D less than 
80, is most probably due to the acoustic excitation present at the end of the third shock 
cell (i.e. x/D = 7). The influence of such excitation on the mean flow seems however to 
be confined to the region x less than 80 D. 

Just as in the previous case, the half width r,/ 2 in the x, z plane grows linearly with x 
for x greater than about 80£>. However, the location of the virtual origin lies downstream 
rather than upstream of the nozzle exit. The crossover point (equal half widths in the two 
central planes) seems to appear much further upstream compared to the cases shown in 
figure 14. The necking down of the jet in the x,z plane can be seen clearly in the Schlieren 
photograph shown in Figure 17b. The fact that y,/ 2 and X|/ 2 form parallel or even divergent 
lines for x greater than 80Z> (see Figures 14 and 15) suggests that true axisymmetry of the 
mean velocity profiles (profiles are identical in both planes) within the region of interest for 
these flows is not obtained, which is in contrast to the observation reported by Trenacosta 
and Sforza 1967, for the case of an incompressible jet where axisymmetry was found for 
x > 150D for an aspect ratio of 10. 

(d) RMS Intensities 
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As explained in the previous section, many uncertainties appear in attempts to deduce 
the velocity fluctuation from a hot-wire voltage in a supersonic flow. It is for this reason 
that the turbulent intensity is defined here as the ratio of the rms of the hot-wire voltage 
to a mean voltage, < e > /E, which allows direct comparison with other work. In 
incompressible flow this is equal to the ratio < u > /U. ; 

Figure 16 shows the downstream variation of the centerline fluctuation intensity 
e = < e > /E 0 for three different pressure ratios. The mean voltage Eo is obtained from 
the relation E 0 = KU 0 , where U 0 is the calculated mean velocity at the center of the 
nozzle exit plane, assuming ideally expanded isentropic flow, and K is a constant obtained 
from the calibration plot of the linearized mean hot-wire voltage versus the mean velocity. 
The value of K is the same for all the data plotted here. Although this approach may 
not correctly represent the precise downstream variation of < u > /£/, the relative trends 
between the various sets of data are expected to be valid. For the underexpanded jet, a 
hot-film probe was used to make measurements for locations of x less than 40D, while 
a hot-wire probe was used for x greater than 40D. For clarity, a line is drawn through 
the data points for the case of 7? = 3.8. Also included in the figure arc the results for a 
high speed subsonic (A/ = 0.8) jet, depicted by the solid line which is an average of the 
data points not shown here. For this case, there is no known acoustic excitation present 
m the flow, thus giving a basis for comparison with the undcrexpanded jet data. As seen 
the magnitude or e increases sharply dose to the jet exit, reaches a maximum value at 
x/D = 10, and then decreases monotonically. Such a behavior is common to most or 
the naturally occurring (non excited) axial turbulence intensity profiles in subsonic jets. 
Generally, a peak in turbulence intensity occurs near the end of the potential core, which 
is a result of fluctuations induced by the interaction of the shear layers. The amplitude 
and position of the peak can be varied significantly by a controlled acoustic excitation of 
the jet near the nozzle exit plane; sec for example Zaman and Hussain 1081. In the self 
similar region of the underexpanded jet, i.c. for x/D > 80, except in the case of It = 3.8, 
the variation of c with x has a behavior similar to that of the subsonic jet. However, the 
magnitude of l at a given location increases with increasing pressure ratio. For the case 
of maximum screech sound radiation, i.e. for R = 3.8, the profile shows a distinct peak 
at x/D = 80, followed by a faster decay in comparison with the two other cases. It was 
found that for 15 < x/D < 50, the magnitude of e decreases like A' -1 / 2 . As discussed 
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below, this type of profile seems to be a result of acoustic excitation in the jet. A weak 
acoustic excitation is also present for the other two pressure ratios (see Figures 2a and 2c), 
and is exhibited here by a small increase in turbulence intensity at x/D — 00, for R = 2.7. 

For an excited plane jet, it has been shown by Hussain and Thomson 1980, that the 
excitation appears to directly afTect the fluctuation intensity and is more pronounced on the 
centerline than in the shear layers. From results for acoustically excited axisymmetric and 
planar jets, it is found that a peak in centerline turbulence intensity occurs in the region of 
interaction of large scale vortical structures, usually appearing near the end of the potential 
core (Hussain and Thomson 1980, Zaman and Hussain 1980). As has already been seen in 
Figure 2b, such structures are present in the region just downstream of the acoustic source. 
Since a peak is observed in the distribution of e at X/D = 80 and in view of the above 
discussion, one would expect to see large-scale vortical structures in a Schlieren photograph 
of the flow field at this location. To explore such a possibility, Schlieren photographs for 
both the x, y and x, z planes, covering the region from the nozzle exit to a downstream 
location of about 110Z?, were taken as shown in Figure 17. Because the available optical 
setup could not view the entire flow field simultaneously, two photographs in each plane 
were taken at separate times using a common reference location and assembled as shown, 
in the figure. Together with the cylindrical wave patteren found in the near sound field of 
the x,y plane (Figure 17a), large scale coherent structures are observed developing along 
the jet boundary in the region of the acoustic source ( also see Figure 2b and 4). Further . 
away from the source no such distinct structures are visible, however, there are large 
random eddies present on the edges of the jet. Figure 17b shows the flow field of the jet 
in the x, r plane, and for x/D greater than 50 a distinguishable large scale structure is 
strongly suggested by the picture. A . dotted line presumably representing such a feature is 
superimposed on the picture. A typical wave length deduced from the picture is found to 
be equal to about 8cm. On comparing a series of randomly selected pictures of this plane, 
there does not appear to be a repeatable pattern in their arrangement. However, the “wave •. 
length” and overall size of these structures match rather well. A more definitive argument , 
concerning the existence of these structures must include quantitative data. With this , 
objective in mind, oscillograph records of hot-wire outputs and their respective frequency , 
spectra were studied. 
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Oscillograph records of e fluctuations along the centerline of the jet, for 
40 < x/D < 140, are shown in Figure 18 for R = 3.8. All the traces have identical vertical 
and horizontal scales. As the probe is traversed downstream, low-frequency oscillations 
dominate the signal, with a recognizable periodicity appearing in the trace at xfD = 80. 
The period can be estimated and is found to be equal to about 1ms (the corresponding 
frequency is 1000Hz), while the signal at x/D = 40, i.e. trace (a), is dominated by 
more random high frequency fluctuations. Following the arguments of Brown and Roshko 
1970, the mean scale of the flow in the x, z plane must increase linearly with (x — xo). 
Consequently, it follows that the mean spacing and the mean size of the eddies or large 
scale structures also increase linearly with (x- xo). The consequence of increasing spacing 
leads to a decreasing frequency, which is consistent with the observations made from the 
oscillograms. 

The passage frequency of the large scale structures at a given location can be obtained 
by studying the frequency spectra of the hot wire signal at that location. Figure 19 shows 
the frequency spectrum of e for various positions along the centerline of the jet for 15 < 
x/D < 140. The horizontal scale is in kilohertz, the range of. which varies for some of the 
plots. For x/D < 40 a strong peak around 17.5kHz is observed and found to be equal to the 
fundamental frequency of the near acoustic field. From a phased locked Schlieren movie of 
the flow field for a 5mm nozzle (Krothapalli et al. 1982), it was found that the frequency 
with which the large scale structures pass any station in the x,y plane was identical to the 
fundamental frequency of the near sound field. Several additional discrete frequencies in 
the range 0-20kIlz are also observed in these plots. With increasing distance downstream, 
the spectral energy shifts toward the low frequencies. At x/D = 80, in addition to other 
low frequency peaks, a distinct peak at 1048Hz is observed, which is close to the frequency 
obtained from trace (c) of Figure 18. This frequency then corresponds to the passage 
frequency of the large scale motions observed in Figure 17b. To further verify this, one 
can calculate the frequency based on the spacing of two adjacent structures in Figure 17b 
and the mean convection velocity. The convection velocity of these structures varies from 
0.5(7* to 0.7 U e . For x/D = 80, we have for the centerline velocity U e =151m/sec; and 
assuming an average convection velocity of 0.6(7*, one obtains a frequency of about 1 1 18llz, 
which is quite close to the measured frequency. It is also interesting to note that the ratio 
of frequencies for the dominant components in the two planes is equal to about 10.7, which 
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happens to be the aspect ratio of the nozzle. 

From thu above it may be suggested that acoustic excitation with sources located at 
x/D < 10 serve to generate the large scale structures (or instability waves) above the line 
scale random background in both the x,y and planes. When the forcing freoueney 

“ " le ' nS,abilUy ba[ld . instability waves (or large scale structures) are triggered 
Winch undergo a process of growth and then breakdown. It has been shown by Crighton 
10.3, that a jet with elliptic cross-section and large eccentricity is virtually stable to all 
temporally growing disturbances: and that in the ease of spatially growing disturbance 
modes representing sideways oscillation parallel to the major axis have a small growth 
rate, while those representing a dapping motion parallel to the minor axis have a large 
growth rate. The disturbances in the above description may be viewed as the vortical 
structures developing within the shear layers of the jet. The growth and break down or 
decay of these disturbances is generally redacted in the amplitude distribution of the rms 
miens, ty along the centerline of the jet. In light of this discussion it is perhaps logical to 
suggest that the amplitude of the disturbances in the *, jt plane increases monolonicnlly, 
reaches a maximum at some distance * less than 15 D, and then rapidly decays, while the 

disturbances in the *, * plane achieve their maximum amplitude at x/D = 80 and then 
decay as shown in Figure 16. 

With regard to the forcing frequency (self excited in this case), one can calculate 
the Strouha! number (5/ = fD/U 0 ) and compare it with the value for any similar flow 
situation to see whether it falls within an instability band. Although the flow near the 
nozzle exit is complicated by the presence of a shock cell structure, we are concerned here ' 
mainly with the flow in the shear layers surrounding the shock cells. This part of the 
flow is mainly compressible and subsonic. Lacking data on the instability characteristics 
of such a flow, it is to a certain extent reasonable to compare the present data with the 
instability characteristics of a plane incompressible air jet. Sato 1060 , experimentally 
investigated the instability of a plane air jot, emerging from a laminar channel flow, in 
the Presence of an imposed transverse acoustic perturbation from a loud speaker placed 
perpendicular to the flow direction. The most unstable jet Strouha! number St =*'/D/U 0 
(D=width of the nozzle) for antisymmetric modes was observed to be equal to O.M. For 
the disturbances in the x, y plane with the observed frequency of 17, 600Hz and using the 
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same definition ( fD/Uo ), the Strouhal number is found to be 0.12, which is within 15% of 
the most unstable Strouhal number of an incompressible plane jet. Similarly, the Strouhal 
number for the disturbances in the plane containing the long dimension of the nozzle, 
i.e. x, z plane, was also calculated, using the observed frequency of 1048Hz and the long 
dimension of the nozzle L =50mm, and found to be again equal to 0.12. Thus, suggesting 
that the most unstable Strouhal numbers in the two planes are the same, except the growth 
rates of the disturbances, as suggested by Crighton 1973, are different in different planes. 
It is also interesting to note that the maximum screech sound radiation occurs when the 
St rouhal number equals to the most unstable Strouhal number of the plane incompressible 
jet. To further verify this, Strouhal numbers for the maximum screech sound radiation 
were calculated using available data in the literature and the results are tabulated, together 
with the present data, in Table 1. Data are also included in the table for axisymmetric 
jets. It appears that the most unstable Strouhal number for rectangular jets (based on 
the criteria of maximum screech sound radiation) appears to be in the range 0.11- 0.14, 
which surprisingly falls within the range of a plane incompressible jet. For axisymmetric 
jets, the Strouhal number (SI = fD/U,D is the diameter of the nozzle) appears to be in 
a range 0.24 - 0.28, which is close to the value 0.3, the so-called “preferred mode” of a jet 
(Zaman & Hussain 1980, Crow & Champagne 1971). 

Profiles of e in the x,y plane at different downstream locations for R = 3.8, are 
shown in Figure 20. In general the profiles shapes are similar to the ones found for two 
dimensional jets, where geometrical similarity exists in the outer flow and a distinct saddle 
shape appears at small values of x/D. For values of x/D greater than 80, the profiles 
tend toward axisymmetry, and seem to approach complete self similarity with increasing 
downstream distance. The profiles for x/D < 80 are much broader in the tj variable in 
comparison with the profiles for x/D > 80. It appears that, when comparing the e profiles 
in Figure 20 with the corresponding mean velocity profiles in Figure 11, the major influence 
of the acoustic excitation is present up to a downstream distance of about 80 D. Profiles of 
e in the x, z plane are shown in Figure 21. For x/D greater 80, it appears that the e profiles 
assume geometric similarity. This observation in conjunction with the self similar mean 
velocity profiles for x/D > 80 suggest that the flow in the z,z plane reaches a self similar 
state. In light of these measurements (both mean and rms profiles) it may be suggested 
that the flow field of a rectangular jet reaches self similar states in both planes at locations 
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downstream of the crossover point. However, the shapes of the self similar profiles are 
cifTerent m the two planes. Furthermore, it appears that complete axisymmetry of the 
jot may not be reached within a reasonable (of engineering interest) distance downstream. 
The profiles of e at different downstream locations for R = 2.7 and 5.4 in both the x,y 
and x, s planes (not shown here), show a similar behavior as discussed above, except for 
the enhanced effects of acoustic excitation. 
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CONCLUSIONS 


/An objective of this work was to study the mean and turbulent velocity fields of an 
underexpanded rectangular jet of moderate aspect ratio. While investigating this problem, 
we encountered two major factors influencing the overall flow field of the jet. One is to 
be expected and is the effect of the density ratio between the jet and the ambient air. 
The second factor, typical of undfcrexpanded jets, is the generation of “screech tones and 
their influence on the structure of the jet. Although discrete tones or “screech" tones are 
found in the near sound field of the jet for all pressure ratios above the critical pressure 
ratio (i.c. R = 1.0), they arc most intense and have greatest effect on the rverall flew field 
only in the range of pressure ratios from 3 to 4. The maximum screech sound radiation 
occurs in this experiment at a pressure ratio of 3.8. When the jet is operating outside this 
range of pressure ratios, the influence of screech tones is less important, as compared lo 
the effects of varying density, on the mean and turbulence velocity Helds of the jet. 

From theSchlieren photographs studied, the following observations are made. For the 
condition of maximum screech sound radiation, a very strong organized cylindrical wave 
pattern is observed, which originates alternately Lorn each side of the jet. The source for 
this wave system is located approximately at the end of the third shock cell and serves as 
an acoustic excitation for the entire jet flow field. Associated with the presence of the wave 
system is an increased spreading rate of the jet in the plane containing the small dimension 
of the nozzle. When the frequency of the acoustic excitation lies within the instability band 
for a plane incompressible jet, it serves to introduce large scale coherent structures in the 
turbulent flow. Such structures are observed here in both planes of the jet. However, 
in the x,y plane they appear just downstream of the acoustic source, i.e. at x/D 10, 
while in the x,z plane they arc observed at downstream locations of x/D greater than 
about GO. The existence of those large scale structures at different downstream locations 
in the two planes was also confirmed by oscillograph records of hot-wire signals and their 
respective frequency spectra. The Stro ibal number for tne.se organized motions was found 
to be equal to 0.12, which is close to the most unstable frequency for the antisymmetric 
modes of the planar incompressible jet. The frequencies in the two planes thus scale with 
the aspect ratio of the nozzle. These observations are supported by the theoretical work 
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of Crighton 107.1, who studied tlu* stability of a jet with an elliptic cross-section. He found 
that in the case of spatially growing disturbance's, modes representing sideways oscillation 
parallel to the major axis, or long dimension of the nozzle, have a small growth rate, while 
those representing a flapping motion parallel to the minor axis, or small dimension of the 
nozzle, have a large growth rate. .... 

The mean and rms intensity profiles, for x/I) greater than 80, in both the central 
x.y and r, r planes, exhibit separate geometrical similarity in their respective planes, thus 
suggesting that complete axisymmetry (i.e. identical profiles in both planes) of the jet 
may not be achieved within a teasonable distance downstream. This is in contrast to out's 
intuitive notion that a rectangular jet shoald reache a complete axisymmctric state in a 
relatively short distance downstream of the nozzle exit. 

In the absence of acoustic excitation with increasing pressure ratio (increase in reser- 
voir pressure), a decrease in the spreading angle of the jet was observed. Similar observa- 
tions have also been made in studios of other supersonic shear layers and is associated with 
the increasing density of the jet. 
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Geometry of Characteristic 
the nozzle dimension D 

Rectangular 0,3cm 

AR=16.7 

Rectangular 0.5cm 

ARsl 0 

Rectangular 0.635cm 

AR=12.75 

Rectangular 0.305cm 

AR=5.oj 

Rectangular 0.18cm 

AR=1 .7 

Axisymmetric 2.54cm 

Axisymmetric 1.27cm 
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Table 1 


Pressure ratio 

Frequency 

Sts 

coiranents 

\ 

velocity U 

Hz 

fd/U 

' V : 


3.8 

17,500 

0.12 

Present 


430 m/s 



results 


3.5 

10,000 

0.12 

Present 


418m/s 



results 


3.67 

7,375 

0.11 

Hanxnitt 


423m/s 





3.67 

18,750 

0.135 

Powell 


423m/s 




; — ' 

3.17 

31,640 

0.14 

Powell 


403m/ s 
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3.67 

4,000 

0.24 


. . ; 

42 3m/ s * 
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9*428 

0.275 

Sherman 



435m/, 
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R= 2 .7 


R= 3-8 


R = 5.4 


Figure 2 


Schlieren pictures of the jet in the i, y plane. Alt = 10.7, l) = 3mm 
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Figure 3. Schlicren pictures of the jet in the x,z plane. Alt = 10.7,/) = 3mm 
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Schlicren pictures of the jet in both x, y and x, z planes 

l/,' = |0, /> = r>mm. 


Figure 4 
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pictures of the jet in the 
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x,y plane, with a reflecting 
AM — 10. P = . r >inm. 


4 \t ; ‘ vj«» k}» 






n * v 1 Ttj ^ lJ 
m V*+,M * An i W 

V d V 1 f ,V if , f kJw* 

H ' V 4 1 

iltol 



f 

h B ■ i 

r| ;/.■•• v;? 








Figure 0. 
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R = 3-7 


Schlicrcn picture of the jet in the r. y plane with a reflecting surface 
placed at 135 to the jet axis. AR = 10, P = 5mm. 
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Figure 13. Mean velocity profiles in the x, z plane, y = 0. R — 3.8; ^ , 

x/D = C>0; 0 , x/D = 80; 0 , x/D = 100; <> , x/D = 120; ,x/D = 
110 . 
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Figure 14 


Growth of a rectangular jet with down stream distance. O • A 


It == 2.7, It = 5.4. 



• • 



• • 








ORIGINAL PAGE is 
OF POOR QUALITY 


Figure 18. Oscillograms of hot wire sig 
H = •'*>: (a) r/D = 10; (b )tjl) 
(c)r/lJ = I •JO: (f) x/[) = i io. 
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